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on the level statistics of a metal grain in a magnetic field 
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A theory is presented for the statistics of the excitation spectrum of a disordered metal grain in 
contact with a superconductor. A magnetic field is applied to fully break time-reversal symmetry in 
the grain. Still, an excitation gap of the order of 5 opens up provided JVP 2 > 1. Here S is the mean 
level spacing in the grain, F the tunnel probability through the contact with the superconductor, and 
N the number of transverse modes in the contact region. This provides a microscopic justification 
for the new random-matrix ensemble of Altland and Zirnbauer. 
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The proximity to a superconductor is known to induce 
a gap in the excitation spectrum of a normal metal. Semi- 
classical theories of this "proximity effect" show that the 
gap closes if time-reversal symmetry (T) is broken (by 
a magnetic field or by magnetic impurities). Recently, 
Altland and Zirnbauer argued that a gap remains in 
the spectrum of a metal grain surrounded by a supercon- 
ductor — even if T is broken completely. (The classical 
mechanics of such a system had previously been studied 
Q.) The gap is small (of the order of the mean level 
spacing in the grain), but it has the fundamental impli- 
cation that the level statistics is no longer described by 
the Gaussian unitary ensemble (GUE) of random-matrix 
theory §. 

The GUE has a probability distribution of energy lev- 
els of the form 

P({E n }) cx H(Ei - Ejf []cxp {-cEl) , (1) 

i<j k 

with some constant c > depending on the mean level 
spacing at the Fermi level (chosen at E = 0). This en- 
semble was first applied to a granular metal by Gorkov 
and Eliashberg and derived from microscopic theory 
by Efetov many years later ||. A single-particle energy 
level E n corresponds to an excitation energy \E n \, that 
is to say, the excitation spectrum is obtained by folding 
the single-particle spectrum along the Fermi level. The 
folded GUE has been studied in Ref. ||. Altland and 
Zirnbauer introduce a different probability distribution, 

P({E n }) cx Y[(E? Eff]jEl cxp (-2cEl) , (2) 

i<j k 

for the (positive) excitation energies of a metal grain in 
contact with a superconductor. (The excitation spectrum 
is discrete for E < A, with A the excitation gap in the 
bulk of the superconductor.) The distribution is re- 
lated to the Laguerre unitary ensemble (LUE) of random- 
matrix theory H by a change of variables. The density of 
states p(E) in this ensemble vanishes quadratically near 
zero energy jl]|7) , 



The gap in the excitation spectrum is of the order of the 
mean level spacing 6. The folded GUE, on the contrary, 
has no gap but a constant p{E) = l/<5 near E = 0. 

In this paper we present the first microscopic theory for 
the effect on the level statistics of the coupling to a super- 
conductor. We consider the case that the conventional 
proximity effect is fully destroyed by a T-breaking mag- 
netic field Q . Assuming non-interacting quasiparticle ex- 
citations, and starting from the well-established GUE for 
the level statistics of an isolated metal grain, we obtain 
a crossover to Altland and Zirnbauer's distribution (||) 
as the coupling to a superconductor is increased. This 
provides a microscopic justification for the "maximum- 
entropy" hypothesis on which Ref. (I]] was based. Such 
a justification is needed because, in contrast to ensem- 
bles in statistical mechanics, there is no physical prin- 
ciple that would require a random-matrix ensemble to 
maximize entropy. Furthermore, because the argument 
of Ref. Q is based on the presence or absence of a cer- 
tain discrete symmetry in the Hamiltonian, it can not 
provide a criterion for how strong the coupling to the su- 
perconductor should be for the new ensemble to apply. 
Our microscopic approach permits us to identify this cri- 
terion, and to compute explicitly how the gap in p(E) 
opens up as the coupling strength is increased. 

We consider the geometry shown in Fig. |l| of a disor- 
dered metal grain (N) , which is connected to a supercon- 
ductor (S) by a point contact or microbridge containing a 
tunnel barrier. Breaking T requires a magnetic field of at 
most a flux quantum through the grain. This field is less 
than the critical field of the superconductor if the size of 
the grain is greater than the superconducting coherence 
length. For simplicity of presentation we consider a real 
order parameter A in S. (We have found that a spatial 
dependence of the superconducting phase, considered in 
Ref. pj , has no effect on the level statistics in the absence 
of T.) We assume zero temperature, so that motion in 
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the grain is totally phase coherent. We seek the distribu- 
tion of the excitation energies £„<A. We first consider 
the density of states p(E). 




FIG. 1. A disordered normal- metal grain (N) coupled to a 
superconductor (S). The black area indicates a tunnel barrier. 

To determine p(E) we adopt the scattering approach of 
Ref. We model the point contact by a normal- metal 
lead supporting N transverse modes at the Fermi level. 
Andreev reflection at the interface scatters electrons into 
holes. This corresponds to the off-diagonal blocks in the 
scattering matrix Sa for Andreev reflection, 



5 A = 











(4) 



where each of the four blocks is an N x N matrix. The 
scattering matrix Sn for the normal-metal grain plus tun- 
nel barrier does not couple electrons and holes, and thus 
has the block diagonal form 



S Q (E) 
SS(-E) 



(5) 



Here So (So) is the scattering matrix for electrons (holes) 
at an energy E from the Fermi level. The N x N scatter- 
ing matrix So can be expressed in terms of the M x M 
Hamiltonian Hq of the isolated grain and an M x N cou- 
pling matrix W 

S (E) = 1 - 2niW f (E -H + iirWW^^W. (6) 

The finite dimension M of Hq is artificial and will be 
taken to infinity later on. 

As demonstrated by Efetov ||, an ensemble of disor- 
dered metal grains in a magnetic field can be described 
by the GUE ffj], 



P(H ) cxexp(-cTr J ff 2N 



(7) 



(Eq. (|l]) follows upon integration over the eigenvectors 
of Ho.) The coefficient c is related to 5 by c = tt 2 /8MS 2 . 
We recall that S is the mean level spacing in the folded 
GUE, which is one half the mean level spacing of Hq. 
The coupling matrix W has the form 



W„ 



(2M5\ 
m = 1,2, 



^) 1/2 (2F« 1 - 1 + 2r « 1 ^ T ^ 



M, n=l,2, 



,N. 



1/2 



(8) 



Here T n is the tunnel probability of mode n through the 
normal lead . For later use we introduce a parameter 
A = 2MS/n and an M x M matrix X = (tt/\)WW t . In 
view of Eq. (^), the matrix X is diagonal with N non-zero 
diagonal elements x n , related to r„ by 



T n 4x n (1 -|- x n ) 



(9) 



The excitation energies E n are the positive roots of the 
equation Det[l — Sa S^{E)\ = 0, which can be rewritten 
as an eigenvalue equation Eq], 



Det(£ -H) = 0, H 



H -XX 
-XX -m 



(10) 



The effective Hamiltonian Ti is the key theoretical inno- 
vation of this work. It should not be confused with the 
Bogoliubov-de Gennes Hamiltonian TLbg , which contains 
the superconducting order parameter in the off-diagonal 
blocks ]l6|. The Hamiltonian TLbg determines the en- 
tire excitation spectrum (both the discrete part below 
A and the continuous part above A), while the effective 
Hamiltonian Ti determines only the low-lying excitations 
E n <C A. As we will see, the spectrum of Ti can be ob- 
tained from a mapping onto a generalization of the well- 
known non-linear cr-model. The Hermitian matrix Ti is 
antisymmetric under the combined operation of charge 
conjugation (C) and time inversion (T): 



H = -C T Ti T C, 



C 



-1 

1 



(11) 



The CT-antisymmetry ensures that the eigenvalues of Ti 
lie symmetrically around E = 0. This discrete symmetry 
(for 7Ybg) was the main point in the maximum-entropy 
argument of Altland and Zirnbauer |lj . 

To compute the spectral statistics on the scale of the 
level spacing, we need a non-perturbative technique. We 
employ the supersymmetric method [p|JlC|] , suit ably mod- 
ified to incorporate the special symmetry (|ll]) of Ti. 
The density of states 



p(E) = ^]imJ-lmF(z) 



is obtained from the generating function 
F(z) = J dj> (exp[±i4> f (E + ie-H + zL) 



( ^ 



Ti = 



Hl 2 






-ft la 



L 



I2M 0-3 
I2MC3 



(12) 

(13a) 
(13b) 
(13c) 



Here is a 4M-component supervector containing 2M 
commuting and 2M anticommuting variables. Half of 
each 2M variables correspond to electron states and half 
to hole states. The charge conjugation operator C inter- 
changes electron and hole variables. The matrices Ti 11 2 
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and I2M 03 are tensor products between a 2M x 2M and 
a 2 x 2 matrix (l p is the p-dimensional unit matrix and 
03 is a Pauli matrix). The appearance of —TL in the CT- 
conjugated block of 71 reflects the CT-antisymmetry ( pd| ) 
of 7i. The measure d(j) is normalized such that F(Q) = 1. 
The brackets (• • •) indicate an average over Hq with dis- 
tribution (Q). 

To evaluate F(z) we perform a series of steps which are 
by now standard in the field |B|,[l(|[l7]] . We first average 
Hq over the GUE, which can be done exactly since it in- 
volves only Gaussian integrals. A term which is quartic in 
4> appears, and we decouple it by a Hubbard-Stratonovich 
transformation. This transformation introduces an addi- 
tional integral over an 8 x 8 supermatrix Q, which we 
evaluate by a saddle-point approximation that becomes 
exact in the limit M — ► 00. We solve the saddle-point 
equation in the limit E — > at fixed N and E/8. As 
in Ref. fl7fl , a manifold of saddle points (determined by 
Q 2 = 1) appears in this limit, while for E ^> 8 only a 
single saddle point remains. 

The matrices Q on the saddle-point manifold have the 
electron-hole block structure 



Qi 
Q 2 



Q 2 = -C T QjC, C = 



CT 3 

la 



(14) 



The 4x4 supermatrix Q\ belongs to the coset space of 
the non- linear a- model in the unitary symmetry class, 
and — Q2 is the CT-conjugate of Q%. The matrix C is the 
charge conjugation operator for the cr-model. The den- 
sity of states is obtained as an integral over the saddle- 
point manifold, 



p{E)=Imi — J dQx Str[Lr(Qi 

xexp[-£ 1 (g 1 )-£ 2 (g 1 )]|, 

1 7T 

= -irA E + le ) Str[r(Qi + Qa) 
Ad 

1 N 

where Str denotes the supertrace and 



Q2 



(15a) 
(15b) 
(15c) 

(16) 



The action C2 can be simplified by expanding it in 
powers of Qi — Q%. This is justified either if r„ <C 1 for 
all n or if N ^> 1. (We therefore exclude the case that 
N and T n are both close to 1.) The first non- vanishing 
term in this expansion is 

-±9A Stv[(Q 1 - Q 2 )% (17) 

N o 2 N nn2 

E ° X 3 = \ ^ 

(1 + X 2 ) 2 ^ 10 - r l2 









fU \ 


7 


03 J 


r = 


V -la J 



£a(Qi) 
5A 



(18) 



The parameter <7a is the Andreev conductance |18|]19| ] of 
the tunnel barrier at the NS interface, which can be much 
smaller than the normal-state conductance g = Ylj=i ^j- 
(Both conductances are in units of 2e 2 /h.) For identical 
tunnel probabilities r,- e T « 1 one has g = NT while 

1 , T ^n 

9A 



I NT 2 . 



Finally, we evaluate the integral ([l5]) using the stan- 
dard decomposition of Qi in terms of angular and radial 
variables ■ The result is 



P(E) 



sm(irE/6) 
~^E 



ds e 



7TE 

T 



As 
9A 



(19) 



Eq. ( |19| ) describes the crossover from p(E) = 1/8 for 
gx < 1 to Altland and Zirnbauer's result (^) for gx 1. 
In Fig. U we have plotted the opening of the gap as the 
coupling to the superconductor is increased. The CT- 
symmetry becomes effective at an energy E for gx ^ E/8. 
For small energies E <C 8 rnin(^/9A, 1) the density of 
states vanishes quadratically, regardless of how weak the 
coupling is. 
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FIG. 2. Density of states for three different values of 
qa = 0.4, 4, 40. The solid curves are the analytical result 
(L9), the data points are from a numerical solution of Eq. 
(LO) [with N = 20, M — 100, and a mode-independent tunnel 
probability Tj = V determined by Eq. (|l8|); some 10 4 random 
matrices Hq in the GUE where generated to compute p(E)]. 
In the inset the analytical result is shown on an expanded 
scale for the same values of qa as in the main plot. The 
dashed line is Eq. (^|), corresponding to the limit <;a —*■ 00. 
The dotted line corresponds to the limit qa — > of a folded 
GUE. 
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As a check on our calculations, we have also computed 
p(E) numerically from the eigenvalue equation (|l0|), by 
generating a large number of random matrices Hq in the 
GUE. The numerical results (data points in Fig. ||) are 
in good agreement with Eq. (fl9|). 

The parameter <?a which governs the opening of the 
excitation gap, does so by enforcing a CT-antisymmctry 
on the non-linear a-model. To see this, consider the 
term (17) in the action, which is proportional to c/a. For 
9a 3> 1 this term constrains Q2 to be close to Qi, and in 
the limit c/a ~> 00 one obtains the CT-antisymmetry 



-C T Q{C = Q 1 . 



(20) 



For (/a <C 1, on the contrary, Q 2 may be quite different 
from Qi, and the CT-antisymmetry is effectively broken. 

We generalized these considerations to level-density 
correlation functions. For this, one has to consider a 
more general source term [replacing the term z L in Eq. 
(|l3|)], and higher-dimensional supervectors (containing 
both advanced and retarded components). After carry- 
ing out the same steps outlined above for the density 
of states, we arrive at a non-linear cr-model with a bro- 
ken CT-antisymmetry. This symmetry is restored for 
<7a - * 00, when the it- model becomes equivalent to that 
associated with the Laguerre unitary ensemble of Ref . [Q . 
This establishes the validity of the distribution (||) in the 
limit of a strong coupling to the superconductor. 

In summary, we have presented a microscopic the- 
ory for the random-matrix ensemble which Altland and 
Zirnbauer obtained from a maximum-entropy hypothe- 
sis. The CT-antisymmetry of the Hamiltonian of non- 
interacting quasiparticles induces an excitation gap even 
if the conventional proximity effect is destroyed by a mag- 
netic field. The Andreev conductance c/a — \ NT 2 of the 
contact between the normal metal and the superconduc- 
tor governs the size of the gap, which becomes of the 
order of the mean level spacing S for c/a 3> 1. An inter- 
esting problem for future research |20[] is the sensitivity 
of the gap to Coulomb interactions between the quasi- 
particles, which break the charge-conjugation invariance 
of the Hamiltonian. 
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